Gluino Air Showers as a Signal of Split Supersymmetry 



Javier G. Gonzalez, Stephen Reucroft, and John Swain 
Department of Physics, Northeastern University, Boston, MA 02115 

It has been proposed recently that, within the framework of split Supersymmetry, long lived 
gluinos generated in astrophysical sources could be detected using the signatures of the air showers 
they produce, thus providing a lower bound for their lifetime and for the scale of SUSY breaking. 
We present the longitudinal profile and lateral spread of G-hadron induced extensive air showers 
and consider the possibility of measuring them with a detector with the characteristics of the Pierre 
Auger Observatory. 
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I. INTRODUCTION 

In the almost structureless fast falling with energy in- 
clusive cosmic ray spectrum, three kinematic features 
have drawn considerable attention for a long time. These 
features, known as the knee, the ankle, and the ultravi- 
olet cutoff, are the only ones in which the spectral in- 
dex shows a sharper variation as a function of energy, 
probably signaling some "new physics" . In particular, if 
cosmic ray sources are at cosmological distances, the cut- 
off is expected at about 10 10 9 GeV, due to the GZK Q 
interactions of the primaries with the microwave back- 
ground radiation. The existence of data beyond the GZK 
cutoff has been puzzling theorists and experimental- 
ists Q , but a clear and widely accepted explanation is yet 
to see the light of day. The proposed resolution for this 
puzzle generally invokes physics from the most favored 
theories beyond the standard model (SM) like string/M 
theory, supersymmetry (SUSY), grand unified theories 
(GUTs), and TeV-scale gravity jj. 

A novel beyond-SM-model proposal to break the GZK 
barrier is to assume that ultrahigh energy cosmic rays 
are not known particles but a new species of particle, 
generally referred to as the uhccron, U 0. The mea- 
ger information we have about super-GZK particles al- 
lows a naive description of the properties of the U. The 
muonic content in the atmospheric cascades suggests 
U's should interact strongly. At the same time, if U's 
are produced at cosmological distances, they must be 
stable, or at least remarkably long lived, with mean- 
lifetime r > 10 6 (my/3 GeV) (d/Gpc)s, where d is the 
distance to the source and mjj, the uhecron's mass. Ad- 
ditionally, since the threshold energy increases linearly 
with mjj, to avoid photo-pion production on the CMB 
m u ^ 1-5 GeV. Within the Minimal Supersymmetric 
(MS) extension of the SM, the allowed range for gluino 
masses is m,g < 3 GeV and 25 GeV < m~ g < 35 GeV. In 
this direction, it was noted in |fj that light Supersymmet- 
ric baryons (made from a light gluino + the usual quarks 
and gluons, mjj < 3 GeV) would produce atmospheric 
cascades very similar to those initiated by protons. 

Recently, Arkani-Hamed and Dimopoulos Q proposed 
an alternative to the MSSM in which the mass spectrum 
of the super-partners is split in two. In this theory, all 
the scalars, except for a fine tuned Higgs, get a mass at a 
high scale of supersymmetry breaking while the fermion's 



masses remain near the electroweak scale protected by 
chiral symmetry. Additionally, all corrections that in- 
volve loops of supersymmetric bosons are suppressed, 
thus removing most of the tunings required to reproduce 
(g — 2) M , B — B mixing and b — > sj 8J. At the same 
time it allows for radiative corrections to the Higgs mass. 
Moreover, very recently it was shown that there exists a 
realization of a "split SUSY" in String Theory HG3. 

An important feature of split SUSY is the long life of 
the gluino due to the high masses of the virtual scalars 
(m s ) that mediate the decay. Indeed, very strong limits 
on heavy isotope abundance require the gluino to decay 
on Gyr time scales, leading to an upper bound for the 
scale of SUSY breaking O(10 13 ) GeV. Additionally, it 
has been pointed out that the detection of gluinos coming 
from astrophysical sources (with trig ~ 500 GeV) leads 
to a lower bound on their proper lifetime of the order of 
100 yr, which translates into a lower bound on the scale 
of SUSY breaking, O(10 n ) GeV HJ. 

In light of this, it is of interest to explore the poten- 
tial of forthcoming cosmic ray observatories to observe 
gluino-induced events. Some signatures of the air showers 
initiated by these long lived gluinos have been presented 
in 0,0]. In this Brief Report, we carry out a more de- 
tailed analysis by generating gluino air showers through 
Monte Carlo simulations and pave the ground for a future 
study on the actual feasibi lity of measuring them at the 
Pierre Auger Observatory The outline is as follows. 
In Sec. II we review the relevant aspects of cosmic ray air 
showers. In Sec. Ill we first carry out Monte Carlo sim- 
ulations of gluino induced showers and then show their 
distinct signatures in the air shower profile and lateral 
spread at ground level. Section IV contains a summary 
of our results. 



II. COSMIC RAY AIR SHOWERS 

When a high energy particle hits the atmosphere it 
generates a roughly conical cascade of secondary parti- 
cles, an air shower. At any given time, the shower can be 
pictured as a bunch of particles, the shower front, travel- 
ing toward the ground at nearly the speed of light. The 
number of particles in the shower multiplies as the front 
traverses the atmosphere, until the particles' energy fall 
below a threshold at which ionization losses dominate 
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over particle creation, after this point the number of par- 
ticles decreases. By the time the front hits the ground, 
its shape is similar to that of a "saucer" with a radius 
that can range from a few meters to a few kilometers. 
The shape of the shower front is actually closer to that 
of a spherical shell, with a curvature of a few kilometers 
for almost vertical showers to more than a hundred kilo- 
meters for inclined ones. The number of particles as a 
function of amount of atmosphere traversed is the "lon- 
gitudinal profile" of the shower. 

The general properties of the longitudinal profile can 
be understood with a simple model and it usually can be 
parameterized by a Gaisser-Hillas function [l4|. 

N e (X) = N„ ( X Xa ) ' (1) 

where N etVaax is the number of particles at shower max- 
imum, Xq is the depth of the first observed interaction, 
I max is the depth at the maximum, and A = 70 g/cm 2 . 
The position of the maximum, X max , depends on the en- 
ergy as well as on the nature of the primary particle. 
With cosmic ray showers however, there are fluctuations, 
associated mainly with the point where the primary first 
interacts, as well as statistical fluctuations in the devel- 
opment of the shower. 

When the shower front reaches the ground it is spread 
over an area of up to a few kilometers. It is then possible 
to study the density of energy deposited (or particle den- 
sities), on the ground as a function of time and position. 
There are a number of parameterizations for such distri- 
butions but they are mostly modified power laws like the 
following Q 
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where r is the distance to the point where the core hits 
the ground, tm is the Moliere radius at two radiation 
lengths above the observation level, a is another empiri- 
cal parameter and r/ is the free parameter that depends 
on the angle. The normalization constant C will depend 
on the energy. 

If the primary particle is a hadron, the first interaction 
will be a hadronic interaction and the number of hadrons, 
mostly pions, will increase with each interaction. Statis- 
tically, in each interaction, about 30% of the energy goes 
into neutral pions that decay, generating an electromag- 
netic cascade. In this way, the energy of the primary is 
split into an electromagnetic part and a "muonic" part, 
that comes from the 7r ± that managed to decay. If the 
incoming primary has a higher energy, the number of in- 
teractions required to lower the energy per particle under 
the threshold at which the pions will most likely decay 
increases, increasing the fraction of the energy that goes 
into the electromagnetic part. As a result, the number 
of muons in a shower scales as E . This in turn im- 
plies that the number of muons for a nucleus of mass 
A relates to the number of muons on a proton shower: 



N A _ jpm N p At this point it ig worth not i n g 

that these numbers are strongly dependent on the par- 
ticular hadronic interaction model used. There are three 
hadronic interaction models commonly used in air shower 
simulations, SIBYLL 0, QGSJET and DPMJET [l7||. 
All of them are extrapolations of models that agree with 
the experimental data but show a different behaviour at 
energies beyond 10 TeV center of mass energy |l8|. 

If the primary particle is a 7-ray the interactions that 
occur are mostly pair production, Bremsstrahlung, ion- 
ization losses and Compton scattering. Also, at energies 
higher than 10 10 GeV, the LPM effect suppresses the 
cross-section for pair production and Bremsstrahlung. 
This results in the shower being more elongated. 



III. GLUINO AIR SHOWERS 

In this section we will study gluino induced showers. 
To carry out this study we use aires, a set of programs 
specifically designed to simulate the extensive air show- 
ers generated by ultra high energy cosmic rays interact- 
ing with the atmosphere. The Aires system is described 
elsewhere and takes into account the relevant 

interactions, including electromagnetic and hadronic in- 
teractions and transport processes. The hadronic model 
used is Sibyll. 

We use a feature of AIRES that allows for the definition 
of special primaries by providing a program that handles 
the first interactions of each primary until the main pro- 
gram (aires) can take over and simulate the rest of the 
shower until it strikes ground. To model the gluino in- 
duced showers we first determine where there will be a 
major interaction and then inject a proton with energy 
equal to the energy deposited by the gluon. We then 
force each proton to have it's first interaction at its cor- 
responding point, giving rise to a hadronic shower that 
is then simulated by the standard [T^, l20j program. 

The gluino containing hadron (hereafter G) cross sec- 
tion is about half the pion-air cross section and the in- 
elasticity is K oc 1 GeV/Mc HG3- 

The gluino mass 

range is not constrained so we could, in principle, study 
it at different scales. According to [V^, the masses acce- 
sible to neutrino detectors are < 170 GeV. In our case we 
try to probe for higher masses, while keeping the fluxes 
within reach. In our simulations we adopt a gluino mass 
of 500 GeV, which yields an inelasticity of 0.002 [HI- 
This small inelasticity is precisely what allows one to 
model a G shower as a series of proton sub-showers sep- 
arated according to the G mean free path, with each 
proton having about 0.002 of the original G energy. 

In our case, this program takes the G-hadron of a given 
energy, mean free path and inelasticity. The nature of the 
particle chosen to be injected at each vertex will deter- 
mine the amount of energy channeled into the hadronic 
shower and the amount of energy going into the electro- 
magnetic shower. For a detector like the surface array 
of the Pierre Auger Observatory, the hadronic part (the 



3 



| Longitudinal development | 















r ° -bt^ 






- <sf it 




r 

: #T 


"... . +++ -v 








: >? 






Tl 


r r T- , T T 1 1 T r , 1 T , , T 1 T . 





500 1000 1500 2000 2500 3000 

Atmospheric depth {g/cm A 2) 



FIG. 1: Longitudinal profile for protons at 10 , 10 18 , 10 19 , 
10 20 eV and a gluino at 5 x 10 18 eV. 

muons) are enhanced over the electromagnetic part (also 
considering a great part of the electromagnetic part has 
died away in flight). This means that, by injecting a pro- 
ton in each vertex we are underestimating the number of 
muons that can be sampled in the ground. 

The proton sub-showers are generated following these 
simple steps, and are repeated until the G-hadron reaches 
ground level: 

• Calculate the point of the next interaction. 

• Decrease the energy of the G-hadron by a factor 
given by the inelasticity. 

• Inject a proton with energy equal to the decrease 
in the G-hadron's energy traveling in the same di- 
rection. 

Once the initial string of protons is generated, the 
aires package takes over simulation of the standard 
physics interactions and transport through the atmo- 
sphere to produce the set of showers. 

In Fig. ^ we present the average longitudinal profile 
for 100 G-hadron induced air showers, along with the 
longitudinal development of a proton with different en- 
ergies. Our shower simulations have been carried out for 
incident zenith angle 75°. It should be clear that the 
development of a G-hadron shower can not be fitted by 
the Gaisser-Hillas function given in Eq. . One of the 
biggest sources of fluctuations in air showers is the first 
interaction point. This means that, in our case, the lon- 
gitudinal development of any particular shower will show 
small fluctuations in the shape, since one of these is com- 
posed of about ten proton showers. It should be noted 
that it might be possible to isolate these events from their 
background even for zenith angles as low as 60°, where 
the atmospheric depth is around 2000 g/cm, still more 
than three times the point where the shower reaches its 
maximum. 

In Fig. |21 we show the arrival time as a function of 
distance to the core along the major symmetry axis of 
the lateral distribution function, again for an average over 
100 showers. A fit to this plot, taking the spherical front 
approximation to determine the radius of curvature of 
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FIG. 2: Time of arrival versus distance to the core, along the 
symmetry axis of the LDF 
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FIG. 3: Average Lateral spread for 10 18 , 10 19 , 10 20 eV proton 
showers and 5 x 10 18 eV gluino showers. 

the shower front, leads to a value around 78 km. In the 
general case of cosmic ray air showers the front is not 
parameterized as a sphere since the curvature is more 
pronounced near the core and flattens out for big radius. 

In the spherical front approximation, the arrival time 
as a function of the distance to the core (r) is, up to 
terms of order greater than r/R 

r r 2 / 2\ 

t(r) = ~-u r ■ u R + — (1 - (u r ■ u R ) J , (3) 

where R is the radius of the front and ur is the unit 
vector pointing in the direction of arrival of the shower 
and u r is the radial unitary vector on the detector plane. 

The lateral distribution function (LDF) for all parti- 
cles also shows a distinct behavior. In Fig. [31 the total 
LDF on the shower plane is plotted, along with the corre- 
sponding ones for protons of different energies. A distinct 
feature of the LDF from G-hadrons is that the slope de- 
pends on the distance to the core, as opposed to proton 
generated ones. It is due to the fact that G-hadron in- 
duced showers are just a superposition of lower energy 
showers with different ages (hence different slopes). The 
younger showers, that are spread over a smaller area, 
show a steeper LDF. The particle densities are on the 
order of the densities from proton showers. These densi- 
ties need to be scaled according to the detector response. 
Also, for very inclined showers, the steepness of the LDF 
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depends strongly on the zenith angle so a search for these 
has to take into account the angular resolution of the ex- 
periment. 

IV. SUMMARY 

Using a simple model to simulate the interaction of 
G-hadrons hitting the top of the atmosphere, we have 
shown that the longitudinal development of a G-hadron 
induced air shower is very distinct from that of protons 
and gamma rays that compose the main background. 
Thanks to the low inelasticity of G-air interactions, the 
G-hadron will produce a sequence of smaller showers of 
almost the same energy, and the value for the cross- 
section, in agreement with previous works, gives a sepa- 
ration between them that makes it impossible to resolve 
them. 

It should then be possible to differentiate between G- 
hadron induced showers and those of the background, 
since they display such a unique profile. In order to assess 
how difficult it will be to actually measure these showers 
we need to consider the characteristics of the detector. 
In the case of fluorescence techniques, a careful study 
of the signal to noise ratio for the intensities associated 



with these showers will be needed in order to correctly 
estimate the correct aperture. 

It was also shown that the lateral distribution of 
particles at ground level is such that it should be mea- 
surable with a ground array. The expected signature, a 
small front radius given by the arrival times, seems not 
to be realized but the LDF does show a varying slope, 
due to the superposition of different age showers. A 
more detailed study also needs to be done in order to 
study the aperture and the discrimination power, since 
these will depend strongly on the characteristics of the 
detector. 

Note added: After this paper was finished new bounds 
on the Mg-SUSY breaking scale plane were reported [2l| . 
Interestingly, a small window for high scale SUSY break- 
ing and Mq = 500 GeV still remains open. 
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